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Within the context of modified gravity and dark energy scenarios of the accelerated universe, we 
study the stability of de Sitter space with respect to inhomogeneous perturbations using a gauge- 
independent formalism. In modified gravity the stability condition is exactly the same that one 
obtains from a homogeneous perturbation analysis, while the stability condition in scalar-tensor 
gravity is more restrictive. 
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The recent discovery that the expansion of the universe 
is accelerated, obtained by studying type la supernovae 
0, and the cosmic microwave background experiments 
showing that the universe has nearly critical density 0, 
call for a theoretical explanation. Two classes of models 
are predominant in the literature: the first class assumes 
that there is a form of dark energy (or quintessence) un- 
clustered at all scales, that accounts for 70% of the en- 
ergy density p of the universe. This dark energy, of un- 
known nature, is necessarily exotic: to generate accel- 
eration in Einstein gravity it must have negative pres- 
sure Pde < —pDs/'i- The best fit to the observational 
data favours an even more exotic dark energy {phantom 
energy or superquintessence) with equation of state pa- 
rameter w = Pde/ Pde < ^1, which is evolving in time 
0. Were a value w < — 1 to be confirmed by observa- 
tions, it can not be explained by general relativity with 
a canonical scalar field 0, which is the most common 
model of dark energy, because of the Einstein-Friedmann 
equation H = —k {P^ + p^) /6 = — K(/)^/2 < 0, which is 
incompatible with < — (equivalent to superacceler- 
ation H > 0). To model an equation of state parameter 
w < — 1, a phantom field Q or a field coupled nonmin- 
imally to the Ricci curvature R (see Refs. 5, 6] for re- 
views) have been considered. These theories can be seen 
as special cases of scalar-tensor gravity, described by the 
action 7] 
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The second class of models does not require the pres- 
ence of exotic dark energy but modifies gravity at large 
scales by introducing non-linear (in R) corrections to 
the Einstein-Hilbert Lagrangian which become dominant 
only at late times (low curvatures) These theo- 

ries often suffer from problems with the post-Newtonian 
limit [ill ll^ or from instabilities j and are not 
yet accepted as completely viable theories, but they are 
nevertheless interesting as the cosmic acceleration that 
we are observing may be the first sign of a departure 
from Einstein's gravity. Furthermore, these models are 
motivated by certain compactifications of M-theory |0| . 

In both dark energy models and modified gravity, de- 
pending on the model adopted, the universe may acceler- 



ate forever or end its existence at a finite time in the fu- 
ture in a Big Rip or sudden future singularitv HI . . 
Such singularities have been classified in Ref. [l5j ; accord- 
ing to this classification, it is known that singularities of 
type I can occur in these models {.l_Sj , but singularities of 
other types are not excluded. 

The fate of the universe depends on the presence and 
size of the attraction basins of attractor solutions in the 
phase space. In many models of both dark energy and 
modified gravity, a de Sitter attractor solution is found 
[21) 1. In this paper we address the issue of the stability 
of de Sitter space in modified gravity and scalar-tensor 
theories. It is straightforward to assess stability with re- 
spect to homogeneous (time-dependent only) perturba- 
tions; however, it is more significant to establish whether 
de Sitter space is also stable with respect to more general 
inhomogeneous (space- and time-dependent) perturba- 
tions. This is a much more difficult task because of the 
gauge-dependence problems associated with this kind of 
perturbations [l^l- In the following we show that for 
modified gravity the stability condition obtained with a 
gauge-independent inhomogeneous perturbation analysis 
reduces to that obtained in the much simpler homoge- 
neous perturbation analysis, whereas this is not the case 
for scalar-tensor gravity. 

We begin from the generalized gravity action 
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which contains scalar-tensor gravity as the case 
/ ((/), R) = ip{(j))R, and modified gravity f{R) when the 
scalar field (j) is absent and fan ^ 0. In the spatially flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 

ds"^ = -dt'^ + a'^it) idx'^ + dy^ + dz^) 

the field equations are 
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where F = df/dR, H = a/ a, and an overdot denotes 
differentiation with respect to the comoving time t. By 
choosing (iJ, </>) as dynamical variables, the equilibrium 
points of the dynamical system Q)-© are de Sitter 
spaces with constant scalar field (iJo,0o)- These solu- 
tions exist subject to the conditions 



where the scalar harmonics Y are the eigenfunctions of 



the eigenvalue problem V^V* Y 



-k^Y. Here h. 



the three-dimensional metric of the FLRW background 
and the operator Vi is the covariant derivative associated 
with hij, while k is an eigenvalue. The vector and tensor 
harmonics Yi and Yij obey 



QHl Fo-fo + 2Vo = 



/o' - 2Vo' = , 



(7) 



Y,^^V,VjY+^Yh,j. (12) 



(8) Bardeen's gauge-invariant potentials 



where Fq = F(0o,i?o), /o = fi^o,Ro), Vo = V (M, 
Vq' = ^ : a prime denotes differentiation with respect 

to 0, and i?o = l^H^ 

Inhomogeneous perturbations of de Sitter space are 
investigated by using the covariant and gauge-invariant 
formalism of Bardeen-EUis-Bruni-Hwang |22l | in the ver- 
sion studied by Hwang 23] for generalized gravity. The 
metric perturbations are defined by 
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are used, with equations similar to eq. (|15|l defining the 
gauge-invariant variables A/, AF, and AR. The first 
order equations obeyed by the gauge-invariant perturba- 
tions are given in Ref. 23] and they simplify considerably 
in the de Sitter background (Hq, (/jq). To first order, they 
are: 
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with 



AR = 6 
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Furthermore, vector perturbations do not have any ef- 
fect to first order in the absence of ordinary matter |2^ 
and de Sitter space is always stable with respect to first 
order tensor perturbations, as can be seen from eq. (|18|l . 
so we only need to worry about scalar perturbations (see 
Ref. [iJl for details). We first consider modified gravity 
theories obtained by setting 0=1 and / = f{R) with 
ffiR ^ in the action (El. The gauge-invariant pertur- 
bations are related by 2^ 
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where a = aQe^°*. By using the fact that 



aoe 

one obtains Ai? 



Ai?, 



evolve according to 
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at late times the term I? ja} can be safely neglected and 
stability is achieved if the coefficient of <i>// in the last 
term of the left hand side of eq. 126|l is non-negative, i.e. 
(using eq. iQ), if 
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The spatial dependence of the inhomogeneous pertur- 
bations is encoded in the eigenvector k of the spherical 
harmonics; the fact that the only term containing k (or 
the physical wave vector kphys ~ k/a) in eq. be- 
comes negligible on a de Sitter background implies that 
the spatial dependence effectively disappears from the 
analysis. Eq. H27|l coincides with the stability condition 
that can be obtained by a straightforward homogeneous 
perturbation analysis of eqs. (0)) and |SJ). Hence, in the 
stability analysis of de Sitter space in modified gravity 
theories one can safely neglect inhomogeneous perturba- 
tions and limit oneself to the much more approachable 
homogeneous perturbations, thus bypassing the gauge- 
dependence problems. However, this conclusion could 
not be drawn a priori but it necessarily relies on the in- 
homogeneous perturbation analysis presented. Further, 
this result has been shown to be true only for the sta- 
bility of de Sitter space and not for different attractor 
solutions that may be present in the phase space. 

Naively, the physical reason for this considerable for- 
mal simplification could be looked for in the fact that. 



during the quasi-exponential expansion of the universe, 
inhomogeneities (and anisotropics 25]) are redshifted 
away; this is not the whole story though, because the 
simplification found for modified gravity does not occur 
in scalar-tensor theories. In fact, in this case, eqs. H18|l . 
(1^ . and (1231, together with 
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if 1 -f 3/^^/(2wo-Fb) 0. The stability condition of de 
Sitter space in scalar-tensor gravity then becomes 
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For general scalar-tensor theories, this condition is more 
restrictive than the corresponding stability condition ob- 
tained from a straightforward homogeneous perturbation 
analysis of eqs. Q)-©, which is 



fo T/" 



< 



However, for scalar-tensor theories of the form 
5 = / d^x 



(31) 



(32) 



with a single coupling function, eqs. (|Sn|l and coin- 
cide li^. The reason for the failure of these equations 
to coincide in the general case can be traced to the right 
hand side of eq. (|16|) . in which perturbations Ai? act as 
a source for the perturbation A0 (roughly speaking, per- 
turbations SH in the Hubble parameter, and their deriva- 
tives, source scalar field perturbations Scj)): such a term is 
absent in the homogeneous perturbation analysis of the 
Klein-Gordon equation Analogously, such a term is 
absent in eq. (|26|l obeyed by the gauge-independent per- 
turbations (or ^a) in modified gravity. This shows 
that although modified gravity is mathematically equiva- 
lent to a scalar-tensor theory 0, 12^ , the corresponding 
physics is not completely equivalent. 
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As an example of application of the stabilit y co ndition 
(|?7jl . let us consider the theory described by |l(llllLl29| 

(33) 

The de Sitter space must satisfy the condition Rq — 
12Hq = fJ,^ and the stability condition (|27|l can never 
be satisfied: this de Sitter space is always unstable. This 
situation can be ameliorated by a quadratic correction - 
in the theory with 

/(i?) =R - ^+aR\ (34) 

the condition for the existence of de Sitter space is again 
Rq = VS /i^ ; by applying the stability condition H27() one 
obtains that de Sitter space is stable if 

and unstable otherwise (in particular for negative a). 
This result agrees with those of Refs. which follow 
from an independent analysis of the effective potential in 
the Einstein conformal frame version of the scalar-tensor 
theory equivalent to the action ^ with f{R) specified 
by eq. (|^ . Furthermore, the theory described by cor- 
rections in both 1/R and R^ has much better chances of 
passing the Solar System tests than the theory based on 
simple 1/R corrections EtI. 

A more complete discussion of the relation between 
stability of modified gravity and of scalar-tensor theo- 
ries, and the application of eqs. (|27|l and (|30() to other 
specific dark energy and modified gravity scenarios will 
be presented elsewhere. 
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